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摘要 :时 间 序 列 土地 利用 时 空 演 变 规律 分 析 是 当前 的 研究 热点 之 一 ,通过 应 用 自 组 织 映射 神经 网 络 施 法 进行 多 时 间 序 列 土 地 利 
用 变化 时 空 一 体 化 表达 与 演变 规律 分 析 ,探索 区 域 土地 利用 变化 模式 。 基 于 北京 而 2005N2907 2009 .2011、2013 年 5 期 土地 利 
用 遥感 分 类 数据 ,构建 自 组 织 映射 神经 网 络 并 利用 其 聚 类 和 降 维 可 视 化 功能 对 室 企 年 份 的 土地 利用 数据 同时 进行 训练 输出 ,发 
现 建设 用 地 耕地、 林地 牧草 地 .园地 的 聚集 模式 ,并 通过 对 输出 神经 元 进行 包 次 聚 类 以 及 土地 利用 变化 轨迹 分 析 ,获得 北京 市 
郊区 5 个 监测 时 相 土 地 利用 变化 的 时 空 演变 特征 。 结 果 揭 示 出 北京 市 郊区 2005 一 2013 年 土地 利用 变化 具有 明显 的 耕地 型 向 
建设 用 地 型 发 展 的 平原 区 演变 特征 ,以 及 向 林地 型 发 展 的 山区 演变 特 很 凤 各 区 的 发 展 具 有 时 间 上 的 顺序 性 ;总 体 上 形成 6 类 
土地 利用 演变 轨迹 。 
关键 词 : 自 组 织 映 射 (SOM) ;土地 利用 变化 ;多 时 间 序 列 ; 时 空 分 析 ; 轨 迹 分 析 
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Abstract: Multiple time series land using spatial-temporal evolution analysis is an important research area. In this study, 
we investigated the spatial-temporal integrated expression of multiple time series land use change. A self-organizing map 
(SOM ) ieural netWork was used to explore regional land use change modes and to analyze what has driven these changes. 
Remote lsensing data for five land use classification data periods (2005, 2007, 2009, 2011, and 2013) for Beijing were 
Used to train the network, and the outputs identified the aggregation modes for building land, farmland, forest land, 
2rasslafd, and gardens by using the clustering, dimension-reducing, and visual functions of the SOM. Then we conducted 
second-step clustering to produce the neuron and build the land use change trajectories that are needed to analyze the 


spatial-temporal features of Beijing suburban land use changes during the five monitoring periods. The results revealed that 
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there were two land use changes in the Beijing suburbs between 2005 and 2013. One was the development of buildings on 
farmland located on the plains and the other was the development of forest land in mountainous areas. Furthermore, 
development in each district had its own time sequences. This meant that we eventually obtained six land use change 


trajectories in total. 


Key Words: self-organizing map (SOM ) ; land use change; multiple time series; spatial-temporal analysis; trajectory 


analysis 


土地 利用 的 时 空 演变 是 指 各 种 土地 利用 类 型 的 空间 分 布 随时 间 的 变化 , 随 着 时 间 序 列 数据 的 增加 ,其 演 
变 分 析 的 难度 和 复杂 度 也 在 增加 ,而 常用 来 分 析 高 维 时 空 数据 的 自 组织 映 射 方法 可 同时 处 理 多 时 间 序 列 数 
据 , 为 多 时 间 序 列 土地 利用 时 空 演 变 规律 与 发 展 模式 分 析 提 供 便 捷 ,为 土地 利用 空间 布局 的 合理 优化 提供 
支持 。 

当前 关于 土地 利用 变化 的 研究 已 取得 了 很 多 成 果 , 杨 国安 等 i 运用 分 形 模 型 通过 对 比 两 个 时 期 北京 市 
各 土地 利用 类 型 的 分 数 维 及 不 稳定 指数 来 研究 北京 市 土地 利用 的 空间 格局 变化 ; 侍 吝 等 ("| 、Long Hualou 
等 5 Fan Q 等 外 均 是 基于 多 期 数据 相 邻 两 期 的 两 两 对 比 的 方式 ,通过 计算 生地 利用 变化 动态 度 .转移 矩阵 等 
方法 来 定量 的 分 析 说 明 土 地 利用 变化 状况 ;Liu H 等 5、 黄 勇 等 1 < 容 菏 芳 等 六 运用 变化 轨迹 分 析 方法 研究 十 
地 利用 空间 格局 的 变化 规律 。 土 地 利用 变化 的 时 空 分 析 包 括 时 间 和 和 空间 两 个 方面 ,当前 的 研究 虽然 多 是 从 时 
空 角度 来 分 析 , 但 对 二 者 的 结合 并 不 是 很 完善 ,针对 空间 研究 其 分 布 模 式 , 针 对 时 间 多 以 时 间 切 片 的 形式 进行 
分 析 。 当 前 对 于 土地 利用 空间 分 布 模式 的 研究 已 比较 成 熟 ,而 对 于 时 间 序 列 数据 的 处 理 多 是 基于 单 期 或 两 期 
数据 之 间 的 比较 ,对 于 超过 两 个 监测 时 相 的 多 时 间 序 列 堆 地 利用 变化 数据 的 时 空 一 体 化 表达 和 对 比分 析 还 有 
竺 于 更 进一步 地 深入 研究 。 同 时 当前 土地 利用 变化 轨迹 多 采用 图 谱 的 表达 方法 , 即 用 轨迹 代码 的 方式 来 表示 
一 种 地 类 向 另 一 种 地 类 的 变化 ,这 种 表达 方式 难 几 直接 对 轨迹 进行 定量 的 分 析 和 可 视 化 表达 。 
自 组织 映 射 (SOM ) 是 一 种 采用 非 监督 式 学 习 对 输入 样本 数据 进行 训练 ,并 将 输入 空间 的 高 维 数据 在 低 
维 (通常 是 二 维 ) 进行 离散 化 表示 的 人 恰 神 经 网 络 。 该 理论 最 早 是 由 芬兰 赫尔辛基 理工 大 学 Kohonen'* 于 
1982 年 提出 的 , 自 组 织 映 射 方法 与 其 他 人 人工 神 经 网 络 方法 的 不 同 之 处 在 于 它 使 用 了 一 个 邻近 函数 来 保持 输 
和信 空间 的 拓扑 性 质 , 从 而 可 以 直接 在 其 输出 面板 上 进行 二 次 聚 类 以 及 轨迹 构建 。SOM 作为 一 种 神经 网 络 聚 
类 算法 常 被 用 来 探索 空间 对 旬 的 聚集 模式 (2 ;另外 SOM 作为 一 种 降 维和 可 视 化 的 工具 也 常 被 用 来 表达 和 
分 析 识 别 潜在 的 模式 2 在 当前 研究 中 ,SOM 的 聚 类 和 可 视 化 的 双重 功能 多 是 结合 起 来 运用 ,并 常 被 各 领 
域 用 来 分 析 处 理 高 维 时 空 数据 ,比如 SOM 已 被 广泛 的 应 用 于 社会 经 济 变 化 “中流 行 病 ' 小 犯罪 (中航 
线 ' 小 气象 .1 等 领域 的 时 空 演 变 和 轨迹 分 析 以 及 可 视 化 研究 中 。 但 到 目前 为 止 ,将 SOM 方法 应 用 到 土地 利 
用 时 空 演变 分 析 中 的 研究 还 比较 缺乏 。 国 外 像 Arribas-Bel 等 .已 开始 尝试 运用 SOM 方法 来 研究 城市 扩张 
及 城市 的 空间 结构 7 而 国内 焦 利 民 等 ' 引 也 开始 尝试 运用 SOM 方法 来 对 中 国 主要 城市 的 扩展 特征 进行 对 比分 
析 , 息 还 证 在 较为 初级 的 阶段 ,多 是 利用 SOM 输出 面板 来 表达 能 反映 城市 扩展 特征 的 各 指标 的 分 布 聚集 状况 
和 差异 ,而 对 于 多 时 间 序 列 数据 的 对 比 以 及 变化 轨迹 的 分 析 和 表达 还 需要 更 进一步 的 深入 研究 。 

本 文 基于 北京 市 2005 .2007 、2009 、2011 、2013 年 5 期 土地 利用 数据 ,运用 SOM 的 神经 网 络 聚 类 功能 将 5 
个 时 间 序 列 数 据 同时 输入 网 络 进行 训练 ,得 到 不 同时 间 序 列 数据 在 时 间 和 空间 上 具有 可 对 比 性 的 输出 聚 类 结 
果 , 并 将 此 结果 在 SOM 输出 面板 以 及 地 理 空 间 中 进行 可 视 化 表达 与 对 比分 析 , 并 借助 SOM 的 拓扑 保持 特性 
在 其 输出 面板 上 构建 变化 轨迹 ,从 而 可 以 直接 定量 的 分 析 和 直观 的 表达 轨迹 的 变化 方向 。 因 此 本 研究 通过 应 
用 SOM 方法 为 多 时 间 序 列 土地 利用 变化 的 时 空 一 体 化 表达 和 对 比分 析 提 供 了 一 种 新 的 尝试 和 探索 视角 ,并 
基于 此 来 深入 分 析 和 揭示 北京 市 郊区 2005 一 2013 年 间 各 区 土地 利用 变化 的 时 空 模式 和 发 展 规律 。 
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1 研究 区 与 研究 数据 


本 文 研究 区 为 北京 市 郊区 的 14 个 区 ,包括 朝阳 区 海淀 区 、 丰 台 区 石景山 区 、 顺 义 区 、 通 州 区 、 大 兴 区 、 房 
山区 .昌平 区 门头沟 区 .怀柔 区 .密云 区 、 延 庆 区 .平谷 区 。 研 究 数据 为 北京 市 2005 .2007、2009 .2011 .2013 年 
5 期 由 SPOT 遥感 影像 解 译 获得 的 土地 利用 分 类 数据 (图 1) ,其 中 除 2005 年 采用 SPOT4 全 色 波 段 影像 和 多 光 

谱 波 段 影像 的 融合 结果 作为 数据 源 ,其 余 均 是 以 SPOT5 多 论 谱 波段 影像 为 数据 源 ， 空间 分 辩 率 为 10m, 分 类 
精度 均 在 82% 以 上 。 研 究 数据 来 自 课题 组 20 年 来 承担 的 北京 市 农村 经 济 研究 中 心 的 长 期 土地 利用 峰 测 项 
目 ,考虑 到 遥感 数据 源 的 一 致 性 以 及 从 2005 年 开始 的 精细 监测 ， | 2005 年 以 后 的 5 个 监测 时 和 进入 
研究 。 由 于 长 期 土地 利用 监测 项 目的 连续 性 一致 性 ,土地 利用 分 类 标准 采用 1984 年 全 国 农业 区 划 委 员 会 制 
定 的 《土地 利用 现状 调查 技术 规程 》 中 的 8 大 类 ,结合 遥感 数据 分 类 特点 最 终 将 地 类 划分 为 了 9 其; 唐 设 用 地 
(将 居民 点 及 工矿 用 地 、 交 通用 地 合并 称 作 建设 用 地 ) 耕地 林地 牧草 地 园地 、 术 域 和 未 利用 土地 。 


国 醒 未 利用 十 地 


图 1 北京 市 2005 一 2013 年 土地 利用 分 类 数据 
Fig.1 Land use classification data of Beijing, 2005—2013 


2 研究 方法 


本 文 的 总 体 技术 流程 如 图 2 所 示 , 主 要 包括 SOM 的 网 络 构建 、 二 次 聚 类 以 及 土地 利用 变化 轨迹 分 析 等 方 
面 的 工作 。 
2.1 SOM 方法 

SOM 的 网 络 是 一 个 两 层 网 络 ,由 输入 层 与 输出 层 ( 又 称 竞争 层 ) 构成 ,其 中 输入 层 用 于 接收 输入 的 训练 样 
本 ,而 输出 层 的 神经 元 一 般 是 按照 二 维 阵列 排列 ,两 层 各 个 神经 元 之 间 实 现 双向 连接 ,SOM 网 络 通过 寻找 最 
优 的 权 值 向 量 即 最 佳 匹 配 神 经 元 来 对 输入 模式 集合 进行 分 类 。SOM 算法 的 步骤 为 :初始 化 各 权 值 向 量 , 即 对 
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输出 层 各 权 值 向 量 赋 小 随机 数 并 进行 归 一 化 ;寻找 输入 


2005、2007、2009、2011、2013 年 


数据 的 获胜 神经 元 ;调整 优胜 邻 域 内 的 权 值 向 量 ;重复 北京 市 土地 利用 遥感 分 类 数据 
寻找 输入 数据 的 获胜 神经 元 及 以 后 的 步 又 ,直到 返 代 终 v 
bE 各 年 份 北京 市 各 区 各 地 类 

止 条 件 被 满足 [e) 用 地 面积 比例 数据 

实验 中 SOM 算法 是 基于 Matlab 编程 实现 , 对 于 了 
SOM 输出 面板 尺寸 参数 的 确定 ,由 于 需要 对 变化 的 轨 Bas al a 
汪 济 行 绘制 视 站 将 合 } 尺寸 设 - v 
OT SOM 输出 面板 的 尺寸 设 a ee 
置 的 足够 大 ,从 而 尽量 使 得 每 一 个 输入 节点 在 输出 的 面 量 的 成 分 平面 图 (获胜 神经 元 ) 
板 空间 中 都 有 只 对 应 于 它 的 单一 获胜 节点 。 本 文 的 输 v 

小 » a > 后 AN 为 | 出 经 元 所 | 迹 b | 
入 数据 为 北京 市 5 个 年 份 14 个 区 的 7 种 地 类 属性 即 输 Re DY ny 
入 数据 共有 5ax14 区 =70 条 ,结合 前 人 研究 经 验 及 多 
次 实验 ,设置 SOM 输出 面板 的 尺寸 为 20x20= 400 , 即 输 土地 利用 空间 土地 利用 时 间 
人 将 聚集 模式 变化 模式 

出 有 400 个 神经 元 ,要 远大 于 输入 的 节点 数 ,SOM 的 网 
络 训练 迭代 次 数 设置 为 10000 次 ,以 保证 训练 结果 的 稳 北京 市 区 级 土地 和 用 思 组 深 开 分 析 结果 
定 ,从 而 获得 各 特征 变量 的 成 分 平面 图 以 及 最 佳 匹 配 神 
经 元 的 位 置 。 图 2 时 总 体 技术 流程 图 
211 SOM 输出 成 分 平面 图 Fig:2 OQverall flow chart of this study 


SOM 输入 数据 为 7 种 地 类 比例 数据 构成 的 7 维 向 
量 ,将 这 7 种 地 类 属性 看 作 是 7 个 特征 变量 ,经 过 网 络 训 练 局 ,输出 神经 元 的 权 值 向 量 即 为 此 7 种 特征 变量 构 
成 。 根 据 各 输出 神经 元 的 位 置 对 每 一 特征 变量 的 值 进行 显 动 ,得 到 各 特征 变量 的 成 分 平面 图 ,可 直观 的 发 现 
各 特征 变量 值 的 聚集 状况 。 

图 3 所 示 为 经 过 网 络 训练 后 在 SOM 输出 面板 止 7 个 地 类 属性 各 自 的 分 布 与 聚集 状况 ,每 张 图 代表 一 种 
地 类 ,颜色 越 深 表示 该 地 类 的 比例 越 高 越 仿 中 ， 可 以 发 现在 SOM 输出 面板 上 建设 用 地 耕地 林地 牧草 地 、 
园地 的 分 布 比较 聚集 ,而 水 域 和 未 利用 土地 的 分 布 比较 杂乱 。 其 中 建设 用 地 比例 的 高 值 聚集 在 SOM 输出 面 
板 的 右上 角 ( 图 3a) ,耕地 比例 的 高 值 聚 集 在 左上 角 ( 图 3b) ,林地 比例 的 高 值 聚集 于 右 下 角 ( 图 3c) ,牧草 地 
比例 高 值 聚集 于 左下 角 ( 图 3d) ,而 园地 比例 的 高 值 聚 集 在 左边 (图 3e) 。 由 此 可 以 看 出 经 过 SOM 网 络 训练 
后 的 输入 数据 在 输入 室 间 串 的 模式 在 输出 空间 中 得 以 识别 和 表达 。 
2.1.2 最 佳 匹配 神经 元 

通过 计算 欧 氏 距离 寻找 距 输 入 数据 最 近 的 输出 神经 元 作为 其 最 佳 匹 配 神经 元 ,根据 输入 数据 在 SOM 输 
出 面板 上 最 佳 匹配 神经 元 的 位 置 ,将 各 年 份 输入 数据 在 输出 面板 上 进行 表示 (图 4) 。 结 合 图 3 中 的 成 分 平面 
图 进行 对 比分 析 , 可 以 从 整体 上 发 现在 本 研究 的 所 有 时 相 内 朝阳 区 ,海淀 区 .石景山 区 ,丰台 区 主要 聚集 在 建 
设 用 地 比例 较 高 的 区 域 ,大 兴 区 .通州 区 顺义 区 一 般 聚集 在 耕地 比例 较 高 的 区 域 ,延庆 区 密云 区 .怀柔 区 、 门 
头 沟 区 等 主要 聚集 在 林地 比例 较 高 的 区 域 ,而 平谷 区 主要 在 园地 比例 相对 较 高 的 区 域 。 
22 本 次 襟 类 

将 输出 的 400 个 神经 元 , 按 其 权 值 向 量 进行 二 次 聚 类 ,二 次 聚 类 采用 k-means 算法 ,经 多 次 试验 确定 聚 类 
为 7 类 时 具有 很 好 的 解释 性 ,并 将 其 分 别 在 SOM 输出 空间 以 及 地 理 空 间 中 进行 可 视 化 。 图 5 中 SOM 输出 面 
板 的 区 域 划分 代表 着 二 次 聚 类 的 结果 ;图 6 中 的 折线 图 表示 每 个 聚 类 结果 的 7 个 属性 值 , 也 即 7 种 地 类 的 面 
积 比例 ,根据 该 比例 确定 各 聚 类 的 土地 利用 结构 类 型 ,其 中 林地 过 渡 型 工 和 开 表 示 其 林地 比例 相对 于 林地 型 
为 低 ,但 仍 占 主要 优势 比例 ,林地 过 渡 型 下 中 建设 用 地 及 耕地 也 占有 不 小 的 比例 ;图 7 是 分 别 将 5 个 年 份 的 聚 
类 结果 在 地 理 空 间 上 进行 可 视 化 ,并 以 多 图 的 形式 来 表现 不 同年 份 之 间 的 变化 。 由 于 是 将 不 同年 份 各 区 域 的 
土地 利用 数据 同时 作为 输入 数据 进行 网 络 训练 ,所 以 在 输出 的 结果 中 ,不 同时 期 ,不同 区 域 的 聚 类 结果 具有 可 
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建设 用 地 比例 /% 
有 61.72—73.81 
图 52.43 一 61.71 
有 44.10 一 52.42 
国 36.98 一 44.09 
国 30.61 一 36.97 
国 23.44 一 30.60 
国 17.80 一 23.43 
国 12.83 一 17.79 
加 8.23 一 12.82 
四 3.64 一 8.22 


耕地 比例 /% 林地 比例 /% 
mm 48.01 一 63.11 ml 68.80 一 83.95 
mm 36.81 一 48.00 用 59.62 一 68.79 
有 27.37 一 36.80 国 52.89 一 59.61 
国 21.94 一 27.36 国 45.92 一 52.88 
国 17.79 一 21.93 国 39.32 一 45.91 
国 14.35 一 17.78 国 32.85 一 39.31 
国 10.85 一 14.34 国 25.56 一 32.84 
国 7.46 一 10.84 国 18.41 一 25.55 
回 4.19 一 7.45 器 12.65 一 18.40 
口 073 一 4.18 口 5.13 一 12.64 
牧草 地 比例 /% 园地 比例 
有 12.90 一 21.69 Bl 12.24 一 18.77 
ml 10.59—12.89 ml 8.94 一 12.23 
国 8.32 一 10.58 ml 5.86 一 8.93 
国 620 一 8.31 ml 4.36 一 5.85 
国 4.59 一 6.19 国 3.60 一 4.35 
国 3.37 一 4.58 国 3.03 一 3.59 
国 2.66 一 3.36 国 2.57 一 3.02 
国 2.00 一 2.65 国 1.95 一 2.56 
器 126 一 1.99 加 1.10 一 1.94 
0.22 一 1.25 口 037 一 1.09 
水 域 比例 /% 未 利用 土地 比例 /% 
ml 3.77 一 8.00 Wl 2.19 3.84 
本 2.23 一 3.76 Em 1.55—2.18 
ml 1.75—2.22 Bm 131 一 1.54 
国 1.49 一 1.74 mm 1.15 一 1.30 
国 1.26 一 1.48 于 0.98—1.14 
国 1.07 一 1.25 上 国 0.81 一 0.97 
国 0.90 一 1.06 站 国 0.61 0.80 
国 071 一 0.89 导 国 0.42 一 0.60 
加 0.51 一 0.70 上 0.20 0.41 
口 028 一 0.50 0.03 一 0.19 


图 3 SOM 输出 成 分 平面 图 
Fig.3 Output component planes of SOM 


对 比 性 。 在 图 5 图 6、 图 7 中 所 表示 的 内 容 均 对 应 着 共同 的 7 个 取 类 ,以 便于 进行 时 空 对 比分 析 。 如 图 6 所 
示 5 各 到 类 的 差异 主要 表现 在 建设 用 地 耕地 林地 方面 的 差异 ,其 次 是 牧草 地 和 园地 ,水 域 和 未 利用 土地 的 贡 
献 较 小 k 
2.3 二 轨迹 分 析 

根据 不 同时 期 各 区 输入 数据 的 最 佳 匹 配 神经 元 在 SOM 输出 面板 中 的 位 置 ,将 各 区 在 不 同时 期 对 应 的 输 
出 点 位 依次 连接 起 来 ,从 而 在 低 维 空间 即 SOM 输出 面板 上 构建 各 区 的 土地 利用 变化 轨迹 ,并 对 变化 轨迹 进行 
聚 类 分 析 ,此 时 将 每 个 区 的 各 个 年 份 的 土地 利用 结构 比例 数据 按照 时 间 顺 序 依 次 排列 组 成 一 个 向 量 ,以 此 来 
表示 该 区 的 变化 轨迹 并 将 其 作为 轨迹 聚 类 的 数据 输入 。 根 据 轨 迹 聚 类 结果 研究 区 域 土地 利用 的 发 展 模式 , 同 
时 将 其 在 SOM 输出 空间 及 地 理 空间 上 进行 可 视 化 。 
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3 北京 市 土地 利用 时 空 演变 分 析 


土地 利用 时 空 演 变 分 析 包 括 空 间 和 时 间 两 个 方面 ， 
土地 利用 空间 聚集 模式 可 以 通过 结合 SOM 输出 面板 中 
各 地 类 属性 的 聚集 状况 和 各 年 份 输入 数据 最 佳 匹配 神 
经 元 的 位 置 来 表达 说 明 ; 而 根据 输出 神经 元 二 次 聚 类 结 
果 可 获取 各 年 份 的 土地 利用 空间 分 布 模式 ,并 通过 多 时 
间 序 列 数据 的 对 比 以 及 变化 轨迹 的 构建 ,分 析 与 可 视 化 
来 说 明 土 地 利用 的 时 空 变化 模式 。 
3.1 土地 利用 空间 聚集 模式 

首先 结合 最 佳 匹 配 神 经 元 的 位 置 以 及 SOM 输出 成 
分 平面 图 中 各 地 类 属性 值 的 分 布 ,可 以 直观 地 发 现 不 同 
年 份 各 区 的 聚集 状况 , 见 2.1.2 节 。 根 据 二 次 聚 类 结 
可 以 定量 的 确定 各 年 份 不 同 区 土地 利用 结构 的 相似 聚 
集 情 况 ,如 图 7。 可 以 发 现在 本 研究 时 段 内 ,北京 市 郊 
区 主要 向 两 个 方向 发 展 ,一 是 东南 平原 地 区 向 建设 用 地 
型 的 发 展 模式 ,表现 为 耕地 型 聚 类 的 减少 及 建设 用 地 / 
耕地 型 聚 类 的 增加 ;一 是 西北 山区 向 林地 型 的 发 展 模 
式 , 表 现 为 林地 过 渡 型 工 聚 类 的 减少 及 林地 型 聚 类 的 
增加 。 
3.2 土地 利用 时 空 变化 模式 
3.2.1 空间 模式 演变 分 析 

(1) 东南 平原 地 区 向 建设 用 地 型 发 展 模式 ,由 图 7 
可 以 看 出 顺义 区 、 通 州 区 、 大 兴 区 由 2005' 年 和 2007 年 
的 耕地 型 转变 为 建设 用 地 /耕地 型 4 其 中 硕 尺 区 转变 较 
早 (2009 年 ) ,通州 区 和 大 兴 区 转变 较 晚 (2013 年 ) 。 该 


EENORE 


a 


图 4 SOM 最 佳 匹配 神经 元 
Fig.4 Best matching unit of SOM 


区 朝阳 区 丰台 区 


ee 
te > 
5000 


六 2005 年 x 2007 年 


4 2009 年 ，。2011 年 、@®@ 2013 年 


Bg 建设 用 地 型 耕地 型 建设 用 地 过 渡 型 
建设 用 地 /耕地 型 [9 林地 过 渡 型 I 
林地 型 EE 中 林地 过 渡 型 [I 


转变 主要 是 耕地 比例 的 减少 以 及 建设 用 地 、 林 地 及 牧草 
地 的 增加 ,这 说 明 顺 义 区 .通州 区 ,大兴 区 具有 作为 城市 
发 展 新 区 的 特征 。 平 原 区 的 朝阳 区 丰台 区 始终 属于 建 
设 用 地 型 ,位 于 山区 过 渡 带 的 海淀 区 、 石 景山 区 始终 属 
于 建设 用 地 过 渡 型 ,该 类 型 的 林地 比例 明显 要 高 。 


图 5 二 次 聚 类 SOM 面板 区 域 划分 
Fig. S Regional division of SOM output plane by second- 


step clustering 


(2 西北 也 区 向 林地 型 发 展 模 式 。 由 图 7 可 以 看 出 房山 区 .昌平 区 .平谷 区 由 2005 年 的 林地 过 渡 型 下 转 
变 为 林地 过 滤 型 和 "延庆 区 密云 区 由 2005 年 的 林地 过 渡 型 I 转变 为 林地 型 ,怀柔 区 .门头沟 区 始终 属于 林地 
型 吐 该 转变 主要 是 林地 用 地 比例 的 逐渐 增加 ,而 耕地 等 其 他 地 类 的 比例 有 所 减少 ,说 明 这 5 个 区 均 具 有 以 林 


地 增长 为 主导 的 土地 利用 变化 趋势 。 


总 的 来 说 ,从 2005 年 到 2013 年 ,平谷 区 .房山 区 .昌平 区 延庆 区 和 密云 区 的 土地 利用 变化 表现 为 从 林地 


过 渡 型 开 到 林地 过 渡 型 [ 再 到 林地 型 的 变化 模式 ,顺义 
建设 用 地 /耕地 型 的 变化 模式 ,而 怀柔 区 .门头沟 区 始终 
型 ,朝阳 区 .丰台 区 始终 属于 建设 用 地 型 。 

3.2.2 土地 利用 变化 轨迹 


x .通州 区 .大兴 区 的 土地 利用 变化 表现 为 从 耕地 型 向 


属于 林地 型 ,海淀 区 .石景山 区 始终 属于 建设 用 地 过 渡 


对 于 14 个 区 所 对 应 的 14 条 轨迹 线 , 通 过 聚 类 来 分 析 轨 迹 方向 的 相似 性 , 本 文采 用 k-means 算法 经 多 次 
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试验 确定 将 14 条 轨迹 聚 为 6 类 时 具有 很 好 的 解释 性 ， 80 
Se es sp 一 建设 用 地 型 
分 别 是 萌 阳 区 .丰台 区 为 一 类 ,海淀 区 .石景山 区 为 。 二 
类 ,顺义 区 .通州 区 .大兴 区 为 一 类 ,昌平 区 .房山 区 为 过 六 oo 一 林地 型 
六 到 型 
类 ,平谷 区 为 一 类 ,怀柔 区 门头沟 区 、 延 庆 区 ,密云 区 为 “关中 So 
一 类 。 各 个 区 的 土地 利用 变化 轨迹 如 图 8 所 示 , 其 中 左 。 朋 8 30 一 林地 过 渡 型 I 
EE 一 建设 用 地 过 渡 型 


图 轨迹 线 的 颜色 与 右 图 各 区 的 颜色 一 一 对 应 , 即 相同 的 
颜色 表示 相同 的 轨迹 聚 类 。 0 
根据 图 8 左 图 的 轨迹 聚 类 图 可 以 发 现 朝阳 区 .丰台 

区 \ 海 诞 区 石景山 区 的 变化 轨迹 均 有 向 SOM 输出 面板 
右上 角 的 延伸 趋势 ,说 明 该 区 域 更 趋向 于 高 建设 用 地 比 ee 

、 时 Fig.6 Percent of different land use type in second-step clustering 
例 的 土地 利用 格局 发 展 , 即 更 趋向 于 中 心 城区 发 展 ; 顺 
义 区 、 通 州 区 、 大 兴 区 趋向 于 向 朝阳 区 ,丰台 区 ,海淀 区 、 


建设 用 地 耕地 ”林地 ”牧草 地 园地 ”水域 未 利用 土地 
用 地 类 型 Land use type 


result 


及 忆 建设 用 地 型 
建设 用 地 /耕地 型 
林地 型 

耕地 型 

记 * ] 林地 过 渡 型 I 
"| 林地 过 湾 型 工 
建设 用 地 过 渡 型 


0 20 40km 


图 7 不 同年 份 的 聚 类 结果 在 地 理 空间 中 的 可 视 化 


Fig.7 Geospatial visualization of clustering result in different years 


景山 区 的 方向 发 展 并 紧 随 其 后 ,表现 为 向 建设 用 地 型 发 展 的 过 渡 阶 段 ;怀柔 区 ,门头沟 区 延庆 区 密云 区 的 
变化 轨迹 均 有 向 SOM 输出 面板 碳 下 角 的 延伸 趋势 ,说 明 这 些 区 域 更 趋向 于 向 林地 比例 高 的 土地 利用 格局 发 
展 , 即 更 倾向 于 向 生态 涵养 功能 的 发 展 ;昌平 区 和 房山 区 的 轨迹 方向 介 于 海淀 区 、 石 景山 区 与 密云 区 、 延 庆 区 
等 之 间 , 而 平谷 区 则 始终 在 园地 用 地 比例 较 高 的 区 域 徘徊 。 结 合 图 5 所 示 二 次 聚 类 结果 的 SOM 面板 ,可 更 加 
直观 地 揭示 各 区 的 土地 利用 变化 模式 。 而 且 轨 迹 线 不 仅 可 以 表现 二 次 聚 类 类 别 间 的 变化 ,还 可 以 表现 某 一 类 
别 中 的 变化 , 像 朝阳 区 .丰台 区 在 建设 用 地 型 中 仍然 向 着 建设 用 地 比例 更 高 的 方向 发 展 。 
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图 8 土地 利用 变化 轨迹 与 可 视 化 


Fig.8 Land use change trajectories and visualization 


4 总 结 与 讨论 


SOM 方法 是 一 种 非常 有 效 的 可 视 化 数据 挖掘 方法 ,可 将 高 维 数据 映射 到 二 维 的 平面 上 ,从 而 更 加 直观 的 
表现 和 识别 潜在 的 模式 。 本 人 研究 运用 SOM 方法 来 探索 分 析 北 京 市 区 级 土地 利用 的 时 空 演变 模式 ,通过 设置 
足够 大 尺寸 的 SOM 输出 面板 ,借助 SOM 的 拓扑 保持 特性 将 多 时 间 序 列 数 据 同 时 输入 SOM 网 络 进行 训练 输 
出 ,从 而 实现 土地 利用 变化 的 时 空 一 体 化 表达 和 对 比分 析 蚊 当前 土地 利用 变化 研究 中 存在 的 对 多 期 数据 处 
理 多 是 基于 单 期 或 相 邻 两 期 两 两 对 比 5 的 问题 提供 了 一 种 新 的 解决 思路 。 

聚 类 与 轨迹 分 析 结 果 与 2006 年 出 台 的 张 京 市 “ 斗 一 五 "功能 区 域 发 展 规划 基本 一 致 ,通过 对 SOM 输出 神 
经 元 进行 二 次 聚 类 及 在 其 输出 面板 上 构建 的 土地 利用 变化 轨迹 ,可 以 发 现 一 些 更 为 精细 的 结果 , 像 朝阳 区 、. 丰 
台 区 海淀 区 石景山 区 均 为 城市 功能 拓展 区 ,但 相 比 与 海淀 区 和 石景山 区 ,朝阳 区 和 丰台 区 更 趋 近 于 中 心 城 
区 的 发 展 ;此 外 ,房山 区 与 昌平 区 虽然 被 规划 为 城市 发 展 新 区 ,当前 仍 以 林地 为 主导 用 地 类 型 , 且 林 地 的 用 地 
比例 有 一 定 的 增长 趋势 ,但 从 轨迹 六 向 上 来 看 也 有 一 定 的 朝 海淀 区 与 石景山 区 的 发 展 趋势 ;而 对 于 顺义 区 、 通 
州 区 和 大 兴 区 ,在 其 区 域 发 展 过 程 中 顺义 区 的 转变 发 生 较 早 。 综 上 可 知 ,北京 市 在 2005 一 2013 年 间 土地 利用 
变化 的 特点 主要 表现 为 建设 用 地 和 林地 的 逐年 增长 以 及 耕地 的 逐年 减少 ,说 明了 城市 在 向 外 扩展 过 程 中 各 地 
类 之 间 此 消 披 长 的 结构 变化 关系 及 建设 用 地 对 农用 地 的 侵占 状况 。 因 此 运用 本 方法 可 以 很 好 的 发 现 区 域 土 
地 利用 变化 的 时 室 模 起 ,同时 也 为 更 多 时 间 序 列 的 土地 利用 变化 数据 的 时 空 一 体 化 表达 和 对 比分 析 提 供 了 一 
种 可 供 选 择 的 更 加 直观 便捷 的 探索 方法 。 

本 研究 的 输入 数据 为 70 条 ,设置 的 输出 神经 元 个 数 为 400 个 ,但 当 研究 区 域 更 大 或 者 对 输入 的 空间 单元 
民 放 划分 更 细 时 ,输入 数据 的 数量 会 大 大 增加 ,此 时 输出 神经 元 的 数量 也 需要 相应 增多 ,也 即 需要 更 大 的 
SeM 输出 面板 的 尺寸 ,这 样 SOM 网 络 训 练 过 程 中 的 计算 量 就 会 大 大 增加 ,这 就 需要 考虑 SOM 的 算法 效率 
问题 。 
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